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ABSTRACT 

We analyze coordinated Hinode XRT and EIS observations of a non-flaring active region to inves- 
tigate the thermal properties of coronal plasma taking advantage of the complementary diagnostics 
provided by the two instruments. In particular we want to explore the presence of hot plasma in 
non-flaring regions. Independent temperature analyses from the XRT multi-filter dataset, and the 
EIS spectra, including the instrument entire wavelength range, provide a cross-check of the differ- 
ent temperature diagnostics techniques applicable to broad-band and spectral data respectively, and 
insights into cross-calibration of the two instruments. 

The emission measure distribution, EAI{T), we derive from the two datasets have similar width and 
peak temperature, but show a systematic shift of the absolute values, the EIS EM{T) being smaller 
than XRT EM{T) by approximately a factor 2. We explore possible causes of this discrepancy, 
and we discuss the influence of the assumptions for the plasma element abundances. Specifically, we 
find that the disagreement between the results from the two instruments is significantly mitigated 
by assuming chemical compositi on closer to the solar photospheric composition rather than the often 
adopted "coronal" composition ()Feldmanlll992D . 

We find that the data do not provide conclusive evidence on the high temperature (logT[K] > 6.5) 
tail of the plasma temperature distribution, however, suggesting its presence to a level in agreement 
with recent findings for other non-fiaring regions. 

Subject headings: Sun: activity; Sun: corona; Sun: UV radiation; Sun: X-rays, gamma rays; Sun: 
abundances; Techniques: spectroscopic 



I. INTRODUCTION 

Understanding how solar and stellar coronae are 
heated to high temperatures is one of the most important 
open issues in astrophysics. Coronal heating is clearly 
related to the strong magnetic fields that fill the atmo- 
spheres of the Sun and solar-like stars, but the mech- 
anism that converts magnetic energy into thermal en- 
ergy remains unknown. Theoretical models based on 
steady heating failed to explain the physical properties 
observed in coronal structures. A promising and widely 
studied framework for understandi ng coron a l heat ing is 
the nanoflarc model proposed by iParkerl (jl972| ). In 
this model, convcctivc motions in the photosphere lead 
to the twisting and braiding of coronal magnetic field 
lines. This topological complexity ultimately leads to 
the formation of current sheets, where the magnetic 
field can be rearranged through the process of mag- 
netic reconnection. This model has been further re- 
fined and adapted to a scenario where nanoflares occur 
in unresolved strands in coronal loops, be low the spa- 
tial resolution of availabl e instrumentatio n (IParkerl 



Cargilllll99l 



Warren et al.l 



Cargill fc K limchuk 1997; KhmchukI 
20031: iParenti et al.ll200Q ). 
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Nanoflare models predict the presence of very hot 
plasma with temperatures in excess of 3 MK in non 
fiaring solar regions; depending on the energy of the 
nanofiarc, temperatures may reach or even exceed 
10 MK. Unambiguous detection of such extreme temper- 
atures in non-flaring solar regions can provide convincing 
evidence for the presence of nanoflares. However, such 
detection is not easy, as the amount of very hot plasma 
produced by nanoflares is expected to be very small. Re- 
cent studies have provided some evide nce of the presence 
of hot plasma in the non- flaring Sun. iMcTiernanI (|2009f ) 
carried out an analysis of the temperature and emission 
measure determined from quiescent plasma during the 
2002-2006 decay phase of solar cycle 23, as derived from 
GOES and RHESSI observations. He found a persis- 
tent faint plasma component with temperatures approxi- 
mately constant during the entire 2002-2006 interval, and 
approximately between 5 and 10 MK. However, GOES 
and RHESSI provided different values of such temper- 
ature, and their results were not necessarily well corre- 
lated; also, this analysis relied on the isothermal plasma 
assumption. 

Other studies tried to identify the hot plasma through 
the X- ray emission obse rved by the X-r ay Telescope 
(XRT : [GoIub et a .1120071 ) onboard Hinode (jKosugi et al.l 
120071 ). iReale et al.l carried out a temperature analysis of 



a non-flari ng active region, fi rst with only XRT multi- 
fihcr data (jReale et al.ll2009bD and then combin ing XRT 



and RHESSI observations (jReale et all l2009aD . Their 
findings point to the presence of small amounts of very 
hot plasmas, with temperatures of ~ 5 — 10 MK, and 
emission measure of the order of few percent of the dom- 
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inant cool component. These characteristics of the emis- 
sion measure distribution are compatible with the predic- 
tions of nanoflarc models. However, they spelled out and 
discussed the main limitations of their study and similar 
analyses. First, XRT is also sensitive to plasma at nor- 
mal active region temperatures (2-3 MK) and thus con- 
tamination from the colder active region plasma is a con- 
siderable obstacle to the detection of the much smaller 
amounts of hot plasma; also, the limited temperature res- 
olution of XRT prevents a detailed study of the cold com- 
ponent. Second, RHESSI sensitivity makes it very hard 
to even detect the quiescent active region plasma. Third, 
instrumental calibrati on is an issue for both instruments. 
iSchmelz et al.l ()2009b[ ) also detected a faint hot temper- 
ature tail to the emission measure distribution of active 
region plasma, and determined its tem perature to be 
around 30 MK. Subsequent analyses by ISchmelz et al.l 
([2009a,) included RHESSI data and, while confirming the 
presence of such hot material, could not reconcile the 
XRT and RHESSI observations using the standard cal- 
ibration of both instruments. A self-consistent solution 
was only found if a series of instrumental parameters and 
the plasma element abundances were adjusted, and the 
temperature of the hot plasma decreased. 

Ample efforts have been devoted to the accurate 
determination of the thermal structuring of coronal 
plasma to derive robust observational constraints to 
the mcchanism(s) of coronal heating. The plasma tem- 
perature distribution of the quiet corona and of active 
regions has been investigated through imaging data 
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2010; Svlwester et al.l l2010( ). Several recent works 



have focused on EUV spectra obtained with the i/m- 
ode Extreme Ultra-v iolet Imaging Spectrometer (EIS; 
iCulhane et al.l [20071 ) which provides good temperature 
diagnostic capability, together with higher spatial 
resolution and temporal cadence than previously avail- 
able (e.g.. IWatanabe et all 120071: iWarren et all 
Patsourakos fc KlimchukI l2009t iBrooks et al.l 
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In the present work, we address the issue of determin- 
ing the temperature distribution of coronal plasma from 
a different perspective: we investigate thermal proper- 
ties of coronal plasma in non-flaring active regions using 
simultaneous Hinode observations with XRT and with 
EIS, which provide complementary diagnostics for the 
X-ray emitting plasma. The multi-filter XRT dataset 
together with EIS spectra, including its entire wave- 
length range, allow to accurately determine the thermal 
structure of the active region plasma, and to explore 
the presence of hot plasma in non-flaring regions. We 
use spectroscopic observations from the Hinode/ElS in- 
strument of a quiescent active region to constrain the 
emission measure distribution of the bulk of the active 
region plasma with the spectral lines observed by EIS 
in t he 171-212A and 245-29lA spectral ran ges (see also 
e.g.. lYoung eralll2007l: IDoschek et al.ll2007l) . Since EIS 
is most sensitive to plasma with temperatures of 0.6- 
2 MK, EIS spectra allow us to accurately determine the 
emission measure distribution of the quiescent active re- 



gion plasma, to evaluate the fraction of the observed 
XRT count rates that it emits, and thus investigate the 
true amount of emission from the nanoflaring plasma. 
Thus, the combination of XRT and EIS observations 
of the same active region allows us to characterize the 
plasma temperature distribution with better detail than 
in previous studies. In fact, while some previous stud- 
ies have made use of data from both imaging and spec- 
troscopic data to constrain the properties of the emit- 
ting plasm a fc.g.. ' Warren et al.ll2010l : ILandi et al.ll2010l : 
lO'Dwver et al. 201^), to our knowledge no previous work 
has carried out a determination of the temperature dis- 
tribution by combining XRT and EIS data, nor a quan- 
titative comparison of the independent analysis from the 
different instruments, as we do here. Independent tem- 
perature analysis from the two datasets provide a cross- 
check of the different temperature diagnostics techniques 
applicable to spectral and broad-band data respectively, 
and insights into cross-calibration of the two instruments. 

The observations are described in Section [2] The data 
analysis and results of the determination of the plasma 
temperature distribution are presented in Section [S] Our 
findings are discussed in Section [4] and summarized in 
Section [S] 

2. OBSERVATIONS 

We observed the non-flaring NOAA active region 10999 
close to disk center, beginning on June 20 2008 around 
23 UT for several hours with both the X-ray Telescope 
and the Extreme Ultraviolet Imaging Spectrometer on- 
board Hinode. The details of the observations are pre- 
sented in Table [TJ 

XRT observed AR 10999 in several filters, with a field 
of view (FOV) of 384" x 384", for about 4 hours start- 
ing June 20 2008 at 23:27 UT, with a cadence of about 
5 minutes in each filter. We analyze observations in the 
following filters: ALpoly, C_poly, Ti_poly, Be-thin, Be- 
med, Al-med. The XRT data were processed with the 
standard routine xrt_prep, available in SolarSoft to re- 
move the CCD dark current, and cosmic-ray hits. Fig- 
ure [1] shows the images of XRT observations in Al-poly, 
Be-thin and Bc-med, integrated over the entire observing 
time (see also Tabled] for details). 

The FOV of the EIS observations analyzed here is 
128" X 128" and was built up by stepping the 1" slit from 
solar west to east over a 3.5 hr period from 23:03 UT 
of June 20th to 02:19 UT of June 21st. The study in- 
cludes full spectra on both the EIS detectors from 171- 
212A and 245-29lA. The exposure time was 90 s and 
the study acronym is HPWOOUULLCCD_RAST. The 
observations were carried out during eclipse season and 
they were not paused during eclipses: the black stripes of 
missing data in Figure [2] correspond to Hinode eclipses. 
The EIS data are processed with the eis_prep routine 
available in SolarSoft to remove the CCD dark current, 
cosmic-ray strikes on the CCD, and take into account 
hot, warm, and dusty pixels. In addition, the radiomet- 
ric calibration is applied to convert the data from photon 
events to physical units. The EIS routine cis_ccd_offset 
is then used to correct for the wavelength dependent 
relative offset of the two CCDs of 1-2 pixels in the X- 
direction, and ~ 18 pixels in the Y-direction. Figure [5] 
shows images obtained from EIS observations by inte- 
grating over narrow wavelength ranges each dominated 
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TABLE 1 

Details of Hinode XR.T and EIS observations analyzed in this paper, and shown in 

Figures [1] and [2] 



XRT EIS 

Al-poly C-poly Ti-poly Bc-thin Be-med Al-mcd 171-212A, 245-29lA 

FOV 384" X 384" 128" x 128" 

START OBS 2008-06-20T23:27:30 2008-06-20T23:03:39 

END OBS 2008-06-21T03:38:50 2008-06-21T02:19:14 

tcxp [s] 4.1 8.2 8.2 23 33 46 90 



Hinode/XRT 20 Jun 2008 23:26 UT 
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Fig. 1. — XRT images of AR 10999, obtained summing all the images taken in a given filter over the ~ 3 hr observation (see Table[l]l: 
Al-poly {left; total integration time tint ~ 131 s), Be-thin (center; ti^t ~ 828 s), and Be-med filter (right; tint ~ 1155 s). 



by a single line with different characteristic temperature 
of formation. Wc also show three small areas of the active 
region which have been selected for the detailed analysis 
of thermal structuring (see next section for details). 

3. ANALYSIS METHODS AND RESULTS 

Inspection of the XRT observations, in all filters, indi- 
cate that the active region is characterized by a modest 
level of variability over a wide range of temperatures (see 
Figure [3]). Therefore, in order to increase S/N, we have 
analyzed the XRT dataset obtained by coaligning the 
images taken in each filter at different times, and then 
summing them up. 

In order to carry out a direct and detailed comparison 
of thermal analysis from XRT and EIS data, we have 
then selected a few subregions. To select these regions 
we have first obtained maps of temperature and emission 
measure over the whole active region, by using estimates 
derived with the so-called combined improved filter ra- 
tio method, devised by Reale et al. (2007). The result- 
ing temperature map is shown in Figure |4l We selected 
three regions, of area approximately 30"x30": EC, which 
in the temperature map homogeneously appears as rela- 
tively cool, and two hotter regions (EH, EHH). The val- 
ues of temperature mean and standard deviation derived 
for these three regions from the CIFR diagnostics are: 
1.5±0.1 MK (EC), 1.9±0.2 MK (EH), and 2.0±0.2 MK 
(EHH). When selecting the boundaries of region EHH 
in the XRT data, we slightly modified the shape (how- 
ever maintaining the area value) in order to avoid the 
contamination spots clearly visible in the temperature 
map as bright blobs immediately to the east and north 



of the selected area. Figure [2] shows the selected regions 
in the EIS Fexv 284A image. For the temperature anal- 
ysis presented in the rest of the paper we integrated the 
XRT and EIS signals over each of these regions, over the 
time intervals when the two instruments were simultane- 
ously observing the selected region. Figure [3] shows that 
the XRT light curves of each region changed little over 
the time period when EIS was observing the same areas. 

Figure [5] shows the EIS full spectra in the three re- 
gions selected for the analysis, with the identification of 
the brightest lines. The EIS exposure time (90 s at each 
location) yields high signal in the strongest lines, mostly 
produced in a temperature range logr[K] 5.8 — 6.3. 
Hotter lines present in the EIS wavelength range have 
generally low intensities, in typical solar non-fiaring con- 
ditions, and arc difficult to detect. We have selected a 
list of lines which are suitable for an accurate derivation 
of the plasma thermal properties. This list (see Table [2]) 
includes strong lines, unblended and with reliable atomic 
data, and we also included hot lines for which however 
we can only derive upper limits from the EIS observa- 
tions of this non-flaring active region. These upper limits 
are nevertheless useful to constrain the high temperature 
component (logT[K] > 6.5). 

In the following we will describe in detail the analy- 
sis methods and results obtained for one of the selected 
regions, EHH. Then we will discuss our findings for all 
three regions. 

3.1. Thermal Structuring from EIS spectra 




Fe XVI \; 





Fig. 2. — EIS images in spectral lines formed at different characteristic plasma temperatures: Hell (256. 32A, log(Tmax[K]) = 4.9), 
Ov (248. 46A, log(rmax[K]) = 5.4), Fevm (185.2lA, log(Tmax[K]) = 5.55), Fexii (195.12A, log(Tmax[K]) = 6.15), Fexv (284.16A, 
log(Tmax[K]) = 6.35), Fexvi (262. 98A, log(rmax[Kl) = 6.4). In the 284A Fexv image we indicate the three regions selected for the 
analysis of the plasma thermal structure (see Figure [4] and @: EC (top left), EH {bottom left), EHH (bottom right). 



For the analysis of EIS spectra and XRT data, we use 
the Package for Interactive Analysis of Line Emission 
(PINTofALE. lKashvap fc Drakd[2000l ) which is available 
on SolarSoft. For the spectral analysis and emission mea- 
sure distribution r econstructio n from EIS spectra we use 
CHIANTI v.6.0.1 (IDere et al.|[l9 97, 2009), with the ion- 
ization balance of Brvans et al.l (l2009i ). Unless explic- 
itly state d otherwise we a ssume "coronal" plasma abun- 
dances of lFeldmanI ()1992| ). 

The observed line intensities provide constraints on the 
plasma temperature distribution, as they depend on the 
plasma emissivities G\{T,n), element abundances Az, 
and differential emission measure distribution DEM (T): 

h^Az j Gx{T,n,)DEM{T)dT (1) 

JT 

where DEM{T) = nl dV/dT [cm-^K-^]. 

Several lines listed in Table [5] are density sensitive, so 
it is necessary to measure the electron density of the 
emitting plasma for each region before calculating the 
contribution functions to be used in the analysis. We 
have used the line intensity ratios of Six, Fexii and 
Fexiv available in the dataset, and the results are listed 
in Tabic [31 The most reliable indicator of the elec- 
tron density is Six. since its lines are no t affected by 
blends and have been found to be reliable bv lYoung et al.l 



(|1998[ ) and iLandi et al.l (|2002[ ). The Fexiv intensity 
ratio is less sensitive to the electron density and pro- 
vides a much coarser estimate; the density values it indi- 
cates are broadly consistent with the Six measurements. 
The Fexii densities, while consistent within the uncer- 
tainties, arc a bit larger than those of Si x, a behav- 
ior that h as already been n oted elsewhere (j Young et al.l 
|2009; Wat anabe et al.ll2009D : this might be due to atomic 
physics proble ms affecting th e Fex ii emissivities (how- 
ever, see also iTripathi et al.l l2010f ). Throughout the 
present work, we used log rie = 9.0 for all regions in the 
calculation of the spectral emissivities from CHIANTI. 

Figure [5] shows the EM loci corresponding to the fluxes 
measured for one of the regions. The curves represent 
the loci of constant flux for a given line, i.e. they are ob- 
tained by taking the ratio of the measured line flux to the 
emissivity function (including the element abundance) of 
that line. The EM loci plot provides a visual indication 
of the temperature range covered by the selected lines 
and passbands, since we are also including EM loci for 
the XRT observations. For the XRT channels the EM 
loci are derived by dividing the measured fluxes by the 
temperature response of that channel, derived assuming 
a given set of element abundances (in the case of Fig. [5] we 
use coronal abundances). Also, if the underlying plasma 
were isothermal the EM loci curves would all cross at a 
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TABLE 2 

Identification and measured fluxes of EIS spectral lines used for 
reconstructing the dem of the three regions selected for the 
analysis and shown in Figure [21 



A [A] 


Ion 


log(Tmax[K]) 


flux 


[erg cm ^ 


s-i] 


Notes'' 




EC 


EH 


EHH 




171.U755 


Fe XIX 


5.9 


o A r\ A o 

5494.52 


loToY.o 


nooo f T 




1 '7/1 t^Q ACl 


Fe X 


6.0 


'7Q in c; v 


12152.6 


VO OO 1 1 






Fe X 


6.0 


4231.15 


7555.12 


4/ob.oo 




ion A noA 


Fe XI 


6.1 


1 niinn a 


1 Q K /I A 

lDv5o4.4 


11523.8 




1 OO 1 (iCiCl 


Fe XI 


6.1 


1 'vo A on 
1 /(54.5U 


unoi /in 
oUzl.4y 


1914.51 




-\ QC or: on 
iOD.OOZU 


Fe XII 


6.2 


an A 1 A~i 
094.141 


1 nnn 

luyy.Do 


Aoa nni 
4zb.UUl 




l5D.<504U 


Fe XII 


6.2 


o 1 on o 1 


5212.59 


onno m 




IOC OQ on 


Fe XI 


6.1 


oolo.by 


QQCQ ni 


a A cic\ i^n 

b4by.bU 




1 no noo c; 


Fe XXIV 


7.3 


Q 1 n non 

oiu.yzy 


r on Q Q 7 


nr\A i^oc; 

v5U4.boo 


u 


1 no Qnon 


Fe XII 


6.2 


A 1 r\'7 c; n 

41U/ .oy 


5858.81 


/I r cn on 




1 ri cr 1 1 on 

lyo.lisU 


Fe XII 


6.2 


1 A nno n 

i4uy2.y 


o 1 om o 

zloUl.z 


16415.0 




or»o r\ A A r\ 

^Uz.U44U 


Fe XIII 


6.2 


1 1 o /I n o 

1124U.Z 


15845.5 


14868.5 




on A cit: AO 


Fe XVII 


6.7 




bU.lool 


A A c; n /I '3 
44.oU4o 


u 


ono i^CiACi 


Caxvi 


6.7 


c: 1 nn a 
O.DiyU4 


-| rr /I O 1 /I 

1 / .4ol4 


DO 'Vl Qn 

OO. / loU 


u 


OAfi nonn 
Z4D.U2UU 


oi VI 


0.0 


'7C\ KQA'7 


ni Q /1 1 c; 
yi.v5410 


OQ ofinn 
ze5.zbUy 




O ACl 1 O ACl 


ol VI 


5.6 


on o /in c; 

zy.z4yo 


Qo cnnQ 


15.6522 




o /I n 1 "Ton 

249.1780 


Ni XVII 


6.5 


on cro^o 




Q 1 O f\ AC\ 

312.949 




250.17U2 


Fe XXII 


7.1 


3.55755 


9.38507 


9.474dU 


u 


253.7880 


Six 


6.1 


379.199 


578.696 


436.514 




258.3710 


Six 


6.1 


2737.30 


4175.65 


3144.01 




261.0440 


Six 


6.1 


1151.82 


1731.62 


1318.26 




262.9760 


Fexvi 


6.4 


173.934 


968.696 


1036.80 




263.7657 


Fe xxiii 


7.2 


14.5618 


12.5289 


9.42194 


u 


264.2310 


Sx 


6.2 


925.718 


975.778 


899.526 




264.7900 


Fexiv 


6.3 


2784.31 


5570.13 


4768.17 




265.0010 


Fexvi 


6.4 


19.6691 


105.484 


99.7671 




270.5220 


Fexiv 


6.3 


1356.52 


3060.05 


2607.03 




272.0060 


Six 


6.1 


1117.27 


1592.69 


1191.82 




272.6390 


Si VII 


5.8 


134.686 


237.159 


99.4883 




275.3540 


Si VII 


5.8 


456.285 


706.011 


317.491 




284.1630 


Fe XV 


6.3 


10317.8 


28644.9 


27965.3 





^ "u" indicates upper limits of fluxes of hot lines which are not actually 
detected in the EIS spectra. 



TABLE 3 

Electron density diagnostic results for each 
region. The density adopted throughout the study 
IS logne = 9.0. Densities are in cm~'^. 



Ion 


Line ratio 


EC 


log Tie 
EH 


EHH 


Six 
Fe XII 
Fexiv 


261.0/253.8 
186.8/195.1 
264.8/270.5 


8 95+0.30 

n 1 5 + 0.15 

9.5±0.4 


9.00±0.25 
9.20±0.20 

g Q + 0.3 


9.10±0.15 

q O+0.3 



single points, therefore this plot also provides some indi- 
cation for the characteristics of temperature distribution 
of the plasma. 

The method of reconstruction of the differential emis- 
sion measure distribution that we adopt for the analysis 
of EIS data runs a Markov-chain Monte Carlo (MCMC) 
algorithm on a set of line fluxes (intensities in broad 
bandpasses can also be used, as for instance in the 
case of XRT datasets; see t j3.2p . and it returns an es- 
tim ate of the DEM that ge nerates the observed fluxes 
(see iKashvap fc Draki 119981 for details on assumptions 
and approximations) . With respect to other methods for 
reconstructing the plasma emission measure distribution 
from a set of fluxes, the MCMC method we adopt has 
the main advantage of providing an estimate of the uncer- 



tainties in the derived DEM. The problem of determining 
the emission measure distribution and its confidence lim- 
its is notorio usly challenging (e .K., ICraig fc Brownlll976l : 
lJudge et all 119971: lJudge(l2010( ). Part of the reason is 
that the emission measure at a given temperature cannot 
be determined independently of temperatures at other 
bins, and therefore the corresponding errors are also cor- 
related. The MCMC method works around this funda- 
mental problem by sampling solutions from the full prob- 
ability distribution of the DEM given the data. It is a 
feature of the MCMC chain that regardless of the con- 
ceptual complexity of the solution space, it can be fully 
explored numerically at a relatively low computational 
cost. Thus, the sampled solutions include the effects of 
statistical noise from the measured data as well as corre- 
lations across temperatures that arise due to overlaps be- 
tween the individual contribution functions of the differ- 
ent lines. From the set of solutions thus obtained, we can 
depict the uncertainty range at temperatures of interest. 
Since it is difficult to display correlations, for purposes 
of clarity, we only show the error ranges computed sepa- 
rately at each temperature bin in the figures. Error bars 
computed for predicted fiuxes, and thus abundances, are 
uncorrelated and thus have the usual meaning. Another 
problem with DEM reconstruction is that the derived 
curves are solutions to a Frcdholm integral equation of 
the first kind, and are thus subject to high-frequency in- 
stability. These instabilities are typically suppressed by 
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Fig. 3. — Lightcurves of XRT observations in all analyzed filters, 
for the three selected regions: EC (top), EH [middle), EHH [bot- 
tom). For each filter the signal is normalized to the intensity in 
the first image. The solid line at the bottom of the plot indicates 
the time interval when EIS was observing the selected region. 



imposing a global smoothness criterion. In the case of the 
MCMC method adopted here however, this restriction is 
relaxed. Smoothing is more physically based, is locally 
variable, and is limited by the width and number of the 
line contribution functions used. This leads to solutions 
that individually have more fluctuations, but on average 
the ensemble of solutions produce an envelope that reli- 
ably determines features that are statistically significant. 
As uncertainties on the measured line fluxes listed in 
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Fig. 4. — Temperature map obtained from the combined im- 
proved filter ratio from the XRT images of AR 10999 shown in 
Figure [T] The three regions selected for the detailed tempera- 
ture analysis are indicated: EC [top left), EH [bottom left), EHH 
[bottom right). When selecting the boundaries of region EHH we 
modified the shape (however maintaining the area value) in order 
to avoid the contamination spots clearly visible in the tempera- 
ture map as bright blobs immediately to the east and north of the 
selected area. 

Table [2] we combine in quadrature the statistical uncer- 
tainties which arc typically very small (of the order of 
few percent), with the uncertainty in the absolute radio- 
metric calibration for EIS which is estimated to be 22% 
(|Lang et all [20061 ) . 

Figure [7] shows the results of the MCMC method to 
derive the emission measure distribution, applied to the 
EIS fluxes measured for the EHH region. In this case, as 
in the rest of the paper, we will plot the emission measure 
distribution EM{T) which is obtained by integrating the 
differential emission measure distribution DEM{T) in 
each temperature bin (here we use a temperature grid 
with constant AlogT = 0.05). The upper panel shows 
the emission measure distribution that reproduces the 
measured fluxes; the associated uncertainties estimated 
as described at the beginning of this section are also plot- 
ted. The bottom panel shows the ratios of predicted to 
observed fluxes for the EIS lines. The comparison of 
measured fluxes with predictions based on the emission 
measure distribution shows that all fluxes are reproduced 
within about 50%, and that the fluxes predicted for the 
hot lines, for which the spectra only provide upper limit, 
are all lower than the upper limits. 

We can investigate the compatibility of this derived 
EM{T) with the XRT observations by computing the 
predicted XRT fluxes and comparing them with the mea- 
sured ones. It is worth noting however that when deriv- 
ing the EM(T) using exclusively the information con- 
tained in the EIS spectra, not all temperature bins are 
well constrained. In particular for these EIS datasets 
which lack measured lines in the high temperature range 
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Fig. 5. — EIS spectra of the three regions selected for the analysis: 
EC (top), EH {middle), EHH {bottom). The identification of the 
spectral lines used for the analysis of the plasma thermal properties 
(listed in Table [2]l are shown in the spectrum in the top panel. 



(logT[K] > 6.5) the hot tail of the EM{T) is poorly 
constrained, but XRT is very sensitive to it, having tem- 
perature responses which peak around logT[K] ~ 7 for 
all filters. Therefore some caution must be applied to 
perform this cross-check. To determine in which temper- 
ature bins the EM{T) is constrained by EIS data we look 
at the contribution of each bin to the line fluxes. Folding 
the EM{T) with the line emissivity we can investigate 
the relative weight of each temperature bin for each line, 
as shown in Figure [S] We then consider as constrained 
the temperature bins where the EMiT) contributes more 
than a threshold (percentage) value to at least one of the 
EIS lines. The choice of the threshold value is somewhat 
arbitrary and we chose a rather conservative 5%. 

Considering the EM{T) in the temperature bins satis- 
fying our selection criterium, we calculate the predicted 
fluxes in the XRT filters (which therefore are strictly 
speaking lower limits to the values predicted by the as- 
sumed EM (T)) which are listed in TablelH together with 
the measured fluxes. The comparison of the two sets in- 
dicates that the EM{T) derived from EIS underpredicts 
the XRT fluxes by a factor ^ 2, the discrepancy be- 
ing slightly larger for the thinner fllters (Al-poly, C-poly, 
Ti-poly) which are more sensitive to the cooler plasma 
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Fig. 6. — EM loci obtained for the measured EIS line fluxes, 
including upper limits (dotted lines) for hot lines in EIS passbands. 
The dashed lines indicate the EM loci for XRT measured fluxes. 
We use different colors for lines of different elements, as indicated 
in the inset. 



TABLE 4 

For region EHH, comparison of XRT measured 

FLUXES WITH THE VALUES PREDICTED USING THE EM{T) 
DERIVED FROM EIS LINE FLUXES (AND SHOWN IN THE 
UPPER PANEL OF FIGURE 0, AND FROM THE FORWARD 
MODEL ING W ITH TWO TEMPERATURE COMPONENTS (SEE 
i|3.2l l. ASSUMING CORONAL ABUNDANCES. 



XRT filter Fobs 
[DN/s] 



prcd.EIS 

[DN/s] 



-^prcdjfm 

[DN/s] 



Al-poly 


122.0 


57.8 (41.8-70.2) 


123.2 


C-poly 


95.8 


44.5 (32.8-53.4) 


101.1 


Ti-poly 


69.1 


31.0 (23.0-36.9) 


65.94 


Be-thin 


11.1 


5.85 (3.69-7.61) 


9.82 


Be-med 


1.59 


1.06 (0.66-1.38) 


1.42 


Al-med 


0.730 


0.483 (0.301-0.628) 


0.673 



^ XRT fluxes predicted using the emission measure dis- 
tribution derived using EIS line fluxes only. The values 
in parentheses represent the range of fluxes predicted by 
the Monte Carlo simulations of the EM(T) which are the 
acceptable emission measure distributions, defining the 
error bars shown in Figure [7] 

^ XRT fluxes predicted using the EM{T) derived from 
the forward modeling with 2T, and using coronal abun- 
dances. 



which is better characterized by the EIS data. 

3.2. Thermal Structuring from XRT data 

We then derive the thermal distribution of the coronal 
plasma exclusively from XRT data, by using two differ- 
ent methods: 

• Forward fitting of distributions of filter ratios of 
several filter pairs, through pixel-by-pixel Monte 
Carlo simulations of the observations, using two 
temperature components (2T); 

• MCMC method, as for the analysis of the EIS spec- 
tra (see ^3T|) . 
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Fig. 7. — Top: Emission measure distribution derived for region 
EHH with the MCMC method (see text) using EIS line fluxes listed 
in Table[2] The error bars provide a measure of the extent to which 
the EM{T) of each temperature bin is constrained by the measured 
fluxes. We plot uncertainties only for the temperature bins for 
which the EM(T) contributes to at least 5% of the flux of at least 
one line (see Figure[8l and text for details). Unconstrained bins are 
plotted with dashed line. Bottom: Comparison of measured EIS 
line fluxes with flux values predicted using the emission measure 
distribution derived with the MCMC method. Different colors are 
used for lines of different elements, in accordance with the colors 
used to show the EM loci (Figure |6]l. 

We use XRT temperat ure res ponses calculated assuming 
"coronal" abundances (jFeldm an 1992) and the new XRT 
calibration which includes a model for the contamination 
layers deposited both on the CCD and on the focal plane 
filters (N.Narukage et al., 2010, Solar Physics, submit- 
ted). For the detailed comparison of XRT and EIS data 
for each of the selected region we only use the subset of 
XRT data taken simultaneously to the EIS observations 
(see Fig.©. 

We have applied the forward fitting method with two 
T components already in previous analysis of XRT data 
aimed at deriving the plasma temperature distribution of 
an active region, in particular to inv estigate the presenc e 
of a high temperature component (jReale et al.l l2009bl ). 
The details of th is method arc thoroughly discussed by 
iReale et al.l ([2009b) . In summary, the pixel fluxes in the 
relevant XRT filters are computed from a basic model 
EM(T) along the line of sight defined by few parame- 
ters (amplitude, central temperature, and width of ei- 
ther single or double top-hat functions). Monte Carlo 
simulations are used to randomize some of the EM(T) 




log T[K] 

Fig. 8. — Each curve represents the fractional contribution of 
each temperature bin to the flux of one of the EIS lines used for 
DEM reconstruction (hot lines for which we only derived upper 
limits are not shown). We use the same colors as in Figure [G] 

parameters and to include Poisson photon noise. This 
procedure is replicated to build whole fake XRT images 
in all filters. The simulated XRT data are then analyzed 
exactly as they were actual data: (1) temperature and 
emission measure maps, i.e., single values for each pixel, 
are obtained for a given filter ratio, and this is done for 
several filter pairs which are sampling the thermal struc- 
ture in different ways; (2) for each considered filter ratio 
an emission measure distribution vs. temperature is built 
by summing up the emission measure values of all pixels 
falling in the same temperature bin. The input EM{T) 
is deemed an adequate representation of the underlying 
plasma emission measure distribution when these curves 
of emission measure distributions vs. temperature are re- 
producing satisfactorily the analogous curves obtained 
from actual data. 

The results of these simulations for region EHH are 
shown in Figure [9l a nd Table IH We no te that, as dis- 
cussed thoroughly in iReale et al.l (j2009b[ ) for the analy- 
sis of XRT observations of another active region, while 
this analysis method is certainly approximate it allows 
us to be sensitive to possible minor EM(T) components, 
e.g., small hot components, and to somewhat constrai n 
them. As for the case presented in IReale et al.l (|2009bf ). 
also here we find that a high temperature plasma com- 
ponent, even if much weaker than the dominating cool 
(logT[K] ^ 6.1 — 6.3), is nevertheless appropriate to re- 
produce the XRT observations. 

As mentioned above, in order to derive the EM{T) 
from XRT data we also use the MCMC method, us- 
ing the flux values integrated in each passband over the 
whole subregion. The errors associated to the XRT 
fluxes are estimated by taking into account the photon 
noise and the calibration uncertainties. Like for the EIS 
fluxes, the statistical uncertainties (we observe bright 
sources and integrate over time to further increase the 
signal-to-noise ratio) are typically very small and the 
errors are largely dominated by the calibration uncer- 
tainties (for an estimate of the latter see Narukage et al. 
2010, Solar Physics, submitted). The resulting errors are 
of the order of ^ 5% for the thin filters, and 15 — 20% 
for thin-Bc and the medium filters. The derived best fit 
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Fig. 9. — Top: Emission measure distribution vs. temperature de- 
rived for region EHH from XRT data. Each curve correspond to a 
filter pair: soft filter ratios (thin solid lines), Al_med/Bc_thin (dot- 
ted line), Bc_mcd/Be_thin (dashed line), Be_med/ Aljned (thick 
solid line). Bottom: Emission measure distribution vs. tempera- 
ture derived for region EHH from XRT data simulated using the 
EM(T) histogram as parent emission measure distribution along 
the line of sight. Lines as in the top panel. 



EM{T) is shown in Figure [TU] (red hne); the associated 
uncertainties are rather large, as should also be expected 
considering the limited constraints that the XRT fluxes 
in six filters can provide to the emission measure distri- 
bution over a large temperature range. 

The differences between the results of the two methods 
are not large, compared to the associated uncertainties. 
The agreement is overall acceptable if one considers that 
the approaches are completely different. The MCMC 
method is model independent but it only uses the total 
region flux values. The forward fitting 2T method makes 
an assumption on the basic EM(T) distribution but it 
keeps the spatial information of the pixel values longer 
in the analysis pipeline, and it is therefore more sensitive 
to small local plasma inhomogeneities (at the scale of 
the single pixel). Moreover this latter method gives large 
weight to the Bc_med/ALmcd filter ratio, and this allows 
to detect the small hot component. 
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Fig. 10. — Comparison of emission measure distribution de- 
rived for regi on E HH independently from EIS lines (with MCMC 
method, see ^3.11 black line), and from XRT, writh two different 
methods: MCMC (red line), and forward modeling with two tem- 
perature components (blue line). For the EM(T) derived using 
the MCMC method we show the estimate of uncertainties for those 
temperature bins which contribute to at least 5% to the flux of at 
least one line or passband (see Figure [8t. 

In Figure [10] we compare the emission measure dis- 
tributions derived independently from XRT and EIS 
data. Considering only the temperature regions where 
the EM(T) are constrained by the data, i.e. the bins 
where the error bars are shown, the different datasets 
yield emission measure distributions that are not too dis- 
similar when the (large) uncertainties are taken into ac- 
count. In particular, they have similar peak temperature 
and width of the cool component. Using the EM{T) de- 
rived from XRT data with the MCMC method to cal- 
culate the expected fluxes for the EIS lines of Table [2l 
wc find that the XRT EM{T) reproduces the measured 
EIS fluxes not too accurately. Specifically, the lines 
with typical temperature of formation logT[K] > 6.3, 
such as Fexiv and Fexv, and the lines formed around 
logr[K] ^ 6.1 — 6.2, such as Fexii and Six, are overpre- 
dicted by a factor ^ 3 (typical range spans approximately 
from a factor 2 to 5); cooler lines (logr[K] < 6.0) are in- 
stead slightly underpredicted (typical range: 0.05 to 
~ 1.5x the measured fluxes). 

3.3. Thermal Structuring from EIS and XRT 

We finally take advantage of all the available infor- 
mation, and derive the plasma EM{T) by applying the 
MCMC method to the measured XRT and EIS fluxes 
together. The EM{T) derived for region EHH by com- 
bining the two instruments is shown in Figure llll com- 
pared to the emission measure distributions derived in- 
dependently from each of the two datasets. When com- 
bining XRT and EIS data, the resulting EM{T) is over- 
all rather similar to the EM{T) derived from EIS data, 
which provide much tighter constraints because of the 
large number of lines, and temperature dependence of 
their emissivity. In the high temperature range however, 
where XRT is more sensitive to the plasma emission and 
EIS provides poor constraints, the XRT-hEIS EM(T) 
departs more from the EIS curve. 

Figure [12] shows how weU the XRT+EIS EM (T) re- 
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Fig. 11. — Comparison of emission measure distribution derived 
for region EHH with diff erent methods: from EIS lines (black line; 
same curve as Fi gur e noj l. from XRT (red line; same curve as curve 
shown in Figure [10} , and using EIS and XRT data together (blue 
curve). 



produces the measured fluxes. While the EIS line fluxes 
are reproduced reasonably well, similarly to the case of 
the EIS EM(T), there is some systematic discrepancy 
between the predicted and the observed fluxes in the 
XRT passbands. Specifically, the three cooler filters, 
Al-poly, C-poly, Ti-poly, are underpredicted by about 
a factor 2. The fluxes in the hotter filters are instead 
well reproduced. We note that for these latter filters, 
which are much less sensitive to cool plasma with re- 
spect to the thin filters, the algorithm of derivation of 
the emission measure distribution can find a satisfactory 
agreement with the observations because of the leeway 
in the determination of the high temperature component 
which is not tightly constrained by the other fluxes, and 
in particular does not contribute significantly to the an- 
alyzed EIS lines. These discrepancies are reminiscent of 
the systematics found when predicting XRT fiuxes us- 
ing the EIS EM{T) (see Table E]), and when using the 
XRT EM{T) to predict the EIS line fluxes (see dis- 
cussion at end of i j3.2p . The cross-check provided by the 
three methods points to a systematic difference with EIS 
data yielding emission measure values a factor ~ 2 lower 
than the values derived from XRT data, as also evident 
from Figure [TTl 

We explored how some of the assumptions we have 
made for our analysis may affect these results. Among 
those, the assumption for the element abundances of the 
coronal plasma has a potentially very significant effect on 
the derivation of the plasma emission measure distribu- 
tion. There is a large body of evidence that element 
abundances in the solar corona do not reflect the so- 
lar ph otospheric composition (e.g.. lMeveilll985t iFeldmanI 
|1992|) . In particular, the chemical fractionation appears 
to be a function of the First Ionization Potential (FIP), 
with low FIP elements, such as Fe and Si, typically 
enhanced in the corona by a factor of a few, whereas 
high FIP elements such as O are thought to have coro- 
nal abundances close to th eir photospheric values (e.g., 
IMevedllQSa IFeldmanI [1991) . 

The intensity of spectral lines is linearly dependent on 
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Fig. 12. — Comparison of measured EIS and XRT fluxes with flux 
values predicted using the emission measure distribution derived 
with the MCMC method using both set of fluxes and shown in 
Figure [TTl We use the same colors as in Figure [6] to indicate lines 
of different elements or XRT passbands. 



the element abundance (assuming that the abundance is 
the same at all temperatures for plasma in the studied re- 
gion; see Eq.[T]). Therefore, if the EM{T) is derived from 
lines all emitted from a single element, e.g. Fe, a change 
in the element abundance will not change the shape of 
the emission measure distribution but merely shift its ab- 
solute value: for instance, a decrease of the abundance 
by a factor 2 will determine an increase of the EM{T) by 
a factor 2. Although our analysis of EIS data includes 
lines emitted by different elements, the large majority of 
the considered lines are emitted by low FIP elements Fe 
and Si. Therefore, if we consider a chemical compos ition 
with a FIP bias different from the IFeldmanI (|1992D set, 
typically adopted for coronal plasma and used for all the 
above analysis, we expect the shape of the EM (T) to not 
change by much qualitatively even its absolute values can 
change significantly. 

For wide band instruments like XRT the effect the 
adopted chemical composition of the plasma on the tem- 
perature response is not as straightforward, as at differ- 
ent temperatures the relative contribution of the emis- 
sion from different elements can change significantly. In 
Figure [T3] we show the effect of the element abundances 
on the temperature response o f the XRT filte rs used 
in this work. We consider the IFeldmanI ()1992[ ) "coro- 
nal" abundances, and a set of abu ndances intermedi- 
ate be tween coronal and photospheric lGrevesse fc Sauvall 
(|1998D . i.e. with a FIP bias of 2 instead of 4 which is typ- 
ically assumed for coronal plasma. For instance for this 
set of abundances Fe, Si and Mg have abundances about 
half the values of IFeldmanI (|1992D . Figure [13] shows that 
in the range logT[K] ^ 6.7 — 7.3 the change in the re- 
sponses is approximately linear with the abundances of 
these low FIP elements which are dominating the emis- 
sion in this high temperature range. At lower temper- 
atures however the XRT filters are affected differently 
depending on the details of their response as a function 
of wavelength and therefore the importance of different 
lines. 

We repeated our analysis using different sets of abun- 
dances, and in particular we derived the EMiT) us- 
ing this set of intermediate abundances described here 

above, and the s et of photospheric abundances of 

iGrevesse fc Sauval (|1998[ ). The results we obtained are 
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Fig. 13. — Ratio of XRT responses, as a function of temperature, 
calculated using intermediate abundances (FI P bias = 2), to the 
responses calculated using coronal abundances llFeldmanlll992l ) . A 
different color is used for each filter used in this work, as indicated 
in the inset. 

sho wn in Figur e [T^ also compared to the EM{T) derived 
with iFeldmanI abundances. The curves have very similar 
shape and are overall compatible within the uncertain- 
ties. However it is clear that, as expected in the light 
of the above discussions, the EM{T) obtained for inter- 
mediate abundances (which are lower than the "coronal" 
abundances) is systematically larger than the one found 
for "coronal" abundances, and the EM(T) derived for 
the photospheric abundances case are even larger. In the 
middle panel of Figure [M] we show the comparison of the 
observed fluxes with the predictions using the EM{T) 
derived for intermediate abundances. Comparing these 
findings with the analogous results of Figure [T^] we note 
that although the overall results are not changed dra- 
matically, the fluxes predicted for the XRT thin filters 
are closer to the observations when using intermediate 
abundances. The increases by about 20-35% are com- 
patible with expectations: the predicted fluxes should 
rise accordingly to the increase of the cool component of 
the EM(T) but be reduced by a factor corresponding to 
the decrease of the XRT response (see Figure [T31) . The 
analogous plot for the case of photospheric abundances 
(bottom panel of Figure [T3]) shows that using these abun- 
dance values the agreement between EIS and XRT im- 
proves even further. 

The values assumed for the element abundances can in 
principle be checked with the EIS spectra. As discussed 
above, the DEM obtained from EIS has been determined 
using almost exclusively lines from low-FIP ions, whose 
coronal abundances are enhanced by a factor of ^ 3 — 4 
over the photospheric values. Therefore, any change in 
the actual abundance of the low-FIP elements from the 
assumed value will result in a systematic shift between 
the observed intensities of the lines from S x to S xiii, 
whose coronal abundance is expected to be close to the 
photospheric one, and their values predicted with the 
DEM. There are ten S lines in the EIS spectra, emitted 
by S X to S XIII, and we have compared their predicted 
and observed intensities to check the corrections to the 
coronal abundances. We found no unambiguous evidence 
of systematic differences in the observed to predicted in- 
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Fig. 14. — Top panel: Comparison of emission measure distribu- 
tions derived with the MCMC method using simultaneously EIS 
and XRT fl uxes, and three different sets of abundances: coronal 
abundances HFeldmanl 1 19921 : black line), intermediate abundances 
(rod line; see text for defi nition of intermediate ab undances), and 
photospheric abundances IIGrevesse fc Sauvallll993 : blue line). For 
intermediate abundances, for which Fe and Si have roughly half 
the abundance as in the coronal set, the EM{T) accordingly in- 
creases; the EM{T) increases even further for photospheric abun- 
dances. Middle and bottom panels: Comparison of measured EIS 
and XRT fluxes with flux values predicted by the emission measure 
distribution derived with the MCMC method (XRT -|- EIS) and: (a) 
interme diate abundances (middle) , or (b) photospheric abundances 
( bottom; IGrevesse fc Sauvailll998l ). 



tensity ratios. Moreover, recent estima tes of the S ab - 
solute abundance have been re vised by iLoddersI ()2003[ ) 
and bv lCaffau fc Ludwi^ (|2007|). who decr e ased t hem by 
a factor 1.5 from the IGrevesse fc Sauvall (I1998D value. 
Since this change is of the same order of the revision 
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of the low-FIP element abundances we propose, we feel 
that the uncertainties in the S predicted line intensities 
arc so large that these ions can not be used to confirm 
the abundance change made necessary by the XRT chan- 
nels in this work. S is in any case a borderline clement 
for FIP effect studies, and abundances of other low FIP 
elements are difficult to determine. While oxygen lines of 
Oiv,v,Vi are present in the EIS wavelength range, they 
are formed at low temperature {\og{T[K]) < 5.5) with 
respect to the hotter coronal EIS lines considered in our 
active region study. In that temperature range the ther- 
mal distribution is not well constrained and therefore the 
abundances cannot be accurately determined. 

While the results shown in Figure [14] seem to suggest 
that the photospheric abundances might be a more ap- 
propriate choice for the active region we studied here, 
we note that the data we have used do not allow us to 
determine the abundances and therefore do not allow us 
to disentangle between the abundance effects and cross- 
calibration issues. We also note that this active region 
was already several days old at the time of the obser- 
vation, and therefore it is reasonable to expect some 
significant effect of chemical fractionation to have oc- 
curred producing departures from the photo spheric com- 
posi tion (see e.g.,IWiding fc Feldmanll2001l : however, see 
also lDel Zannall2003D . 

We have presented the complete analysis carried out 
for one of the selected regions, providing insights into 
the thermal distribution of the plasma, limitations of the 
methods, and cross-calibration of the Hinode XRT and 
EIS. Before discussing in details these findings in [2] we 
briefiy present and discuss the analogous results found for 
the other two selected regions, EC and EH, comparing 
them with the results obtained for region EHH. 

Carrying out the same analysis described above m3.ll 
13.21 13. 3|) on the other two regions, we find results quali- 
tatively similar to our findings for region EHH discussed 
above, in particular in terms of the comparison between 
different analysis methods and systematic discrepancies 
between XRT and EIS. 

In Figure 1151 wc show the emission measure distribu- 
tions obtained for the three regions, using the MCMC 
algorithm applied to XRT and EIS data together {top 
panel), or forward fitting the XRT observations using 
two temperature components (bottom panel). We note 
that despite the forward fitting 2T model adopts consid- 
erable simplifying assumptions for the underlying emis- 
sion measure distribution the results we obtain with this 
method are qualitatively in good agreement with the 
MCMC approach which does not impose constraints on 
the shape of the DEM. On the other hand, we note 
that the forward fitting 2T model has the advantage of 
also taking into account the spatial information, while 
the MCMC method only uses the information contained 
in the integrated fluxes in each passband. The two 
hotter regions, EH and EHH, have very similar under- 
lying EM{T): their cool components have analogous 
peak temperature (logr[K] ~ 6.2 — 6.3) and amounts of 
plasma, whereas the hot component of EHH is of weight 
comparable to the high temperature component of EH 
but it appears shifted towards slightly higher tempera- 
tures. For region EC the cool component is character- 
ized by a peak temperature slightly cooler than the other 
two regions, and in particular its EM{T) falls off faster 



on the high temperature side of the cool peak; the hot 
component of EC is much weaker than the other two 
components, if at all present. 



E 10 




5.5 6.0 6.5 

log T[K] 



7.0 




5.5 



6.0 



6.5 



7.0 



Fig. 15. — Emission measure distributions for the three selected 
regions: EC (red line), EH (blue line), EHH (black line). Top: 
EM(T) curves derived with MCMC method, using EIS and XRT 
fluxes. Bottom: EM{T) deriv ed f rom XRT data through the for- 
ward fitting 2T method (see i|3.2l for details), using photospheric 
abundances. 



4. DISCUSSION 

We have presented the results of a detailed analysis 
of Hinode XRT and EIS observations of a non-flaring 
active region to diagnose the temperature distribution of 
the coronal plasma. In this work we carried out indepen- 
dent temperature analysis from measurements of either 
instruments, and use the derived thermal distributions 
to make prediction for the fluxes observed with the other 
instrument and then compare them with actual measure- 
ments. Then we combine the information from both 
instruments and compare the results with the findings 
based on XRT or EIS data only. This approach allowed 
us to explore the limitations of the data and analysis 
methods in providing constraints to the plasma temper- 
ature distribution, and investigate the cross-calibration 
of the two Hinode instruments. With respect to previous 
works using both XRT and EIS data to stu dy the ther- 
mal properties of the coronal plasma (e.g., iLandi et al.l 
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|201C1[) , in this work we provide a quantitative cross-check 
between the diagnostics of the two instruments. 

When using the two datasets separately the derived 
EM{T) curves have overall similar width and peak tem- 
perature. However, the emission measure derived from 
EIS is systematically smaller, by a factor ~ 2 — 3, than 
the emission measure derived from XRT fluxes, for each 
of the three studied regions. We note that even if the 
uncertainties in the derived EM{T) are rather large the 
systematic discrepancy appears to be significant. We find 
that the extent of the discrepancy between XRT and EIS 
depends in part on the assumptions made for the chem- 
ical composition of the X-ray emitting plasma: it can 
be mitigated by assuming abundances i ntermediate be- 
tween the typical "coronal" a bundances (|.F eldman"1992D 
and the photospheric values iGrevesse fc Sauval (1998), 
and it almost disappears when using photospheric abun- 
dances. This interesting finding encourages further in- 
vestigations to search for similar evidence in other active 
regions, in an active region at different epochs, and in 
different kinds of coronal structures (quiet sun, bright 
points, coronal holes). While our study provides no ro- 
bust constraints on the element abundances we stress the 
importance of the assumed values in the analysis of Hin- 
ode observations. The detailed EM{T) reconstruction 
and the instrument cross-calibration are significantly af- 
fected by the assumed abundances: EIS spectra allow 
the determination of "abundanc e independent" EM(T ) 
e.g. using exclusively Fe lines (IWatanabe et al.ll2007f ). 
but its absolute values depend on the abundances; the 
XRT passbands include significant contribution of sev- 
eral elements, both high-FIP and low-FIP, and therefore 
the dependence of the EM{T) on abundances is more 
complicated and temperature dependent. 

For the analysis of XRT data we use two different 
methods for deriving the emission measure distribution. 
We find that the Markov-chain Monte Carlo method 
yields results in good qualitative agreement with the for- 
ward fitting 2T method which makes rather simplifying 
assumptions on the emission measure distribution but 
has the advantage of making use also of the spatial in- 
formation. This finding lends further confidence to the 
results obtained previously by applying this method to 
the study of XRT observations of another active region 
(|Reale et al.ll2009bf ). 

We find that the EIS spectra allow an accurate deter- 
mination of the cool (logT[K] < 6.5), most prominent, 
component, but when used alone the EIS data are unable 
to constrain the hotter emission due to the lack of strong 
hot lines in the spectra. The combination with XRT 
data provide much tighter constraints on the emission 
measure distribution on a wider temperature range. Our 
analysis shows that the studied active region is character- 
ized by bulk plasma temperatures of 2 MK, which are 
rather low for active region plasma. The EM{T) of the 
three selected regions are very similar for logT[K] < 6.2, 
but for the two hotter regions, EH and EHH, the cool 
component is broader and with a peak shifted towards 
higher T. The high temperature tail of the emission mea- 
sure distribution is not strongly constrained by the data. 
However, the EA4 (T) derived at least for the two hotter 
regions, EH and EHH, suggests the presence of a hotter 
component (logT[K] > 6.5) about two orders of magni- 
tude weaker than the dominant cool component, in good 



agreement with the findings of iReale et al.l (|2009bD for 
another, hotter, active region. 

5. CONCLUSIONS 

We analyzed Hinode XRT multi-filter data and EIS 
spectral observations of a non-fiaring active region to 
study the temperature distribution of coronal plasma, 
EM{T), and to carry out a detailed investigation of the 
cross-calibration of the two instruments. We selected 
three subareas of the active region, and for each of them 
we derived the emission measure distribution EM{T) 
by: (1) using EIS measured line fluxes; (2) using XRT 
fluxes in six of the instrument's fllters; (3) combining the 
datasets for the two instruments. 

We find a good consistency in the qualitative character- 
istics - peak temperature and width of dominant temper- 
ature component- of the EM{T) derived with different 
methods. However, the emission measure distributions 
derived by the two instruments XRT and EIS indicate 
a systematic discrepancy between the two instruments, 
with EIS data yielding EM (T) consistently smaller, by 
about a factor 2, than the EM{T) compatible with XRT 
data. We discuss the possible origin of the disagreement 
and find that the assumptions for the element abun- 
dances significantly infiuence the plasma temperature di- 
agnostics. In particular we find that a chemical compo- 
sition intermediate between the us ually adopted coro- 
nal abundances by iFeldmanI (|1992| ) and the solar pho- 
tospheric abundances improves the comparison between 
the results obtained with the two instruments. When 
adopting photospheric abundances the discrepancy be- 
tween EIS and XRT decreases further. However we note 
that it seems unlikely that the observed plasma, in an ac- 
tive region which is not newly emerged, has photospheric 
composition. Furthermore, the used data do not allow a 
definite determination of the abundances, and therefore 
do not allow us to robustly assess of the cross-calibration 
of the instruments. 

One of the main aims of this work was to exploit 
the complementary diagnostics for the X-ray emitting 
plasma provided by XRT and EIS, to investigate the 
presence of hot plasma in non-flaring regions and test 
nanoflare heating models. We find that the derived 
EM{T) are characterized by an expected dominant cool 
component (typically logr[K] ^ 6.3), and a much weaker 
amount of plasma at higher temperature. While the 
amount of hot plasma is in general in agreement with 
recent findings for other non-fiaring active regions, and 
it is compatible with expectations from nanoflare mod- 
els, we find that within the uncertainties these results are 
not conclusive. 
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